Effective spin systems in coupled micro-cavities 
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We show that atoms trapped in micro-cavities that interact via exchange of virtual photons can 
model an anisotropic Heisenberg spin-1/2 chain in an external magnetic field. All parameters of 
the effective Hamiltonian can individually be tuned via external lasers. Since the occupation of 
excited atomic levels and photonic states are strongly suppressed, the effective model is robust 
against decoherence mechanisms, has a long lifetime and its implementation is feasible with current 
experimental technology. The model provides a feasible way to create cluster states in these devices. 
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Introduction: Interacting two level systems, either 
termed spins or qubits, are of central importance in quan- 
tum information and condensed matter physics. Two or 
higher dimensional magnetic compounds arc believed to 
host some of most interesting condensed matter phenom- 
ena, such as frustration and high T c superconductivity 
[lj|, which are not yet fully understood. In quantum in- 
formation, lattices of interacting spins can be employed 
to generate highly entangled states, such as cluster states 
which are the required resource for one way quantum 
computation While it is a prerequisite for quantum 
information processing, the ability to address individual 
spins in an experimental device can also be very helpful to 
obtain deeper and more detailed insight into condensed 
matter physics. 

In magnetic compounds where spin lattices appear nat- 
urally the addressability of individual spins is unfortu- 
nately extremely hard to achieve because the spatial sep- 
aration between neighboring spins is very small and the 
timescales of interesting processes can be very short. 

Here we show that effective spin lattices Q can be gen- 
erated with individual atoms in micro-cavities that are 
coupled to each other via the exchange of virtual photons. 
Due to the size and separation of the micro-cavities, in- 
dividual lattice sites can be addressed with optical lasers, 
whereas the cavities can be arranged arbitrarily allowing 
for various lattice geometries. The two spin polarizations 
| |) and | J,) are thereby represented by two longlived 
atomic levels of a A level-structure (c.f. figures Q] and 
[2]). Together with external lasers, the cavity mode that 
couples to these atoms can induce Raman transitions be- 
tween these two long-lived levels. Due to a detuning be- 
tween laser and cavity mode, these transitions can only 
create virtual photons in the cavity mode which medi- 
ate an interaction with another atom in a neighboring 
cavity. With appropriately chosen detunings, both the 
excited atomic levels and photon states have vanishing 
occupation and can be eliminated from the description. 
As a result, the dynamics is confined to only two states 
per atom, the long-lived levels, and can be described by 



a spin-1/2 Hamiltonian. Due to the small occupation 
of photon states and excited atomic levels, spontaneous 
emission and cavity decay are strongly suppressed. All 
these results arc verified by detailed numerics. A rcaliza- 




FIG. 1: Level structure, driving lasers and relevant cou- 
plings to the cavity mode to generate effective a x a x - and 
CT^o-^-couplings for one atom. The cavity mode couples with 
strengths g a and g^ to transitions \a) <-» \e) and |6) <-> |e) 
respectively. One laser with frequency ui a couples to transi- 
tion \a) <-> |e) with Rabi frequency Q a and another laser with 
frequency u)b to \a) <-> |e) with Qb- The dominant 2-photon 
processes are indicated in faint gray arrows. 

tion of the scheme thus requires cavities that operate in 
a strong coupling regime with a high cooperativity factor 
and an atom photon coupling that exceeds cavity decay. 
Such regimes have now been achieved in several devices 
0, 0, @| , making a realization of the presented scheme 
feasible with current technology. We begin by showing 
how to engineer effective a x a x , a v a y and a z a z interac- 
tions as well as the effective magnetic field a z and then 
explain how to generate the full anisotropic Heisenberg 
model. We verify the validity of all approximations by 
comparison with the exact dynamics of the full atom- 
cavity model and also apply the model to the generation 
of cluster states. Finally we discuss the feasibility of our 
model for realistic experimental parameters. 
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XX and YY interactions: Wc consider an array of 
cavities that are coupled via exchange of photons with 
one 3-level atom in each cavity (figure []} . Two long lived 
levels, | a) and \b), represent the two spin states. The 
cavity mode couples to the transitions \a) <-> |e) and 
| b) *-> \e), where |e) is the excited state of the atom. 
Furthermore, two driving lasers couple to the transitions 
| a) «-> |e) respectively \b) <-» |e). The Hamiltonian of the 



atoms reads Ha = SjL 



w e e 



c-ji t uj a b\bj)(bj\, where 
the index j counts the cavities, uj e is the energy of the 
excited level and uj ab the energy of level \b). The energy 
of level | a) is set to zero and we use h = 1. The Hamil- 
tonian that describes the photons in the cavity modes 



TT t i 7 V"*JV 

is H c = wcL F i fl tt j + J cL,= 



a]a j+1 + aja] +1 



where creates a photon in cavity j, u>c is the en- 
ergy of the photons and Jc the tunneling rate of pho- 
tons between neighboring cavities Q. For convenience 
we assume periodic boundary conditions, where He 
can be diagonalized via the Fourier transform a k 



N 



k = IT ' -f ^ 1 ^ f (^ odd ) t0 



VlV ^J=l ^ "J ' iV ' 2 - " — 2 

give iJc = J2k w fe°l a fc wi th ^fc = w c + 2Jccos(fc). Fi- 
nally the interaction between the atoms and the pho- 
tons as well as the driving by the lasers are described 
by H AC = Ef =1 [(l^e-^* + 9*^) |e,) ( aj \ + h.c] + 
+ [a «-> 6]. Here g a and f?b are the couplings of the 
respective transitions to the cavity mode, Q a is the 
Rabi frequency of one laser with frequency uj a and fib 
the Rabi frequency of a second laser with frequency 
LOb 01- The complete Hamiltonian is then given by 
H = Ha + He + U ac ■ 

We now switch to an interaction picture with respect 



to H = H A + H c - 5t Ef=i \bi){bj 
uj a b — (u) a — ujb)/2, and adiabatically eliminate the ex 



where <5i 



cited atom levels |ej) and the photons 0. We consider 
terms up to 2nd order in the effective Hamiltonian and 
drop fast oscillating terms. For this approach the de- 
tunings A Q E W e - LJ a , A 6 EE LU e - LU b - (u} ab - (Si), 

Sa = w e - uj k and 5% = lu c - uj k - (u ab - 5 X ) have 
to be large compared to the couplings fl a ,fl b ,g a and 
g b , i.e. |A a |,|A b |, I^U^I » |n |, |fi 6 |, \g a \, \g b \ (for all 
k). Furthermore, the parameters must be such that the 
dominant Raman transitions between levels a and b are 
those that involve one laser photon and one cavity pho- 
ton each (c.f. figure [T]). To avoid excitations of real 
photons via these transitions, wc furthermore require 



5 k n » 



^a9b 

2A a 



2A,, 



(for all fc). 



Hence whenever the atom emits or absorbs a virtual 
photon into or from the cavity mode, it does a transition 
from level \a) to \b) or vice versa. If one atom emits 
a virtual photon in such a process that is absorbed by 
a neighboring atom, which then also does a transition 
between |o) to \b), an effective spin-spin interaction has 
happened. Dropping irrelevant constants, the resulting 



effective Hamiltonian reads 

JV 



3=1 



where crj 



\bj}{bj\ - \a 3 )(aj\ and er+ = \bj)(a,\. The 
parameters B, J\ and J 2 arc given to second order by 
[l0(. If J£ = J 2 , this Hamiltonian reduces to the XY 
model, 



H xy = ^ j Ba z j +J x a x j a 



j'+i 



(1) 



with J x = (Ji + J 2 )/2 and J y = (J x - J 2 )/2. 

For fl a = ±(A a g a / A b g b )fl b with fl a and fif, real, the 
interaction is either purely a x a x (+) or purely a y a v (— ) 
and the Hamiltonian (jTJ) becomes the Ising model in a 
transverse field, whereas the isotropic XY model (J x = 
J y ) [11] is obtained for either fl a — > or fib — > 0. The 



effective magnetic field i? in turn can, independently of 
Ja, and J y , be tuned to assume any value between \B\ 3> 
l^xli \ Jy\ an d |-B| -C I J x \i \ Jy\ by varying <5i. Thus we will 
be able to drive the system through a quantum phase 
transition. Now we proceed to show how to engineer 
effective ZZ interactions. 

ZZ interactions: To obtain an effective o~ z o~ z interac- 
tion, we again use the same atomic level configuration 
but now only one laser with frequency to mediates atom- 
atom coupling via virtual photons. A second laser with 
frequency v is used to tune the effective magnetic field via 
a Stark shift. The atoms together with their couplings to 
cavity mode and lasers are shown in figure O The Hamil- 
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FIG. 2: Level structure, driving lasers and relevant couplings 
to the cavity mode to generate effective cr z (j z -couplings for 
one atom. The cavity mode couples with strengths g a and gb 
to transitions \a) <-> \e) and \b) <-> |e) respectively. Two lasers 
with frequencies u and v couple with Rabi frequencies Q a 
respectively A a to transition \a) <-> |e) and Q b respectively A b 
to 1 6) <-» |e). The dominant 2- photon processes are indicated 
in faint gray arrows. 

tonians Ha of the atoms and He of the cavity modes thus 
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have the same form as above, whereas Hac now reads: 



Hac = 



N 



' 2 

and f2;, 



+ 9aaj) \ej){aj 



h.c. 



. =1 [(" 2 

+ [a <-> b] . Here, O a and O^, are the Rabi frequencies 
of the driving laser with frequency cu on transitions 
| a) — * \e) and \b) — ► |e), whereas A a and A b are the 
Rabi frequencies of the driving laser with frequency v 
on transitions \a) — ► |e) and \b) — ► |e). 

We switch to an interaction picture with respect to 
Hq = Ha + He and adiabatically eliminate the excited 
atom levels \ej) and the photons |9| . Again, the detunings 
A a = uj e - uj, A b = uj e - uj - uj ab , A a = uj e - v, A b = 
u) e -v - uj ab , 6% = u) e - uj k and 6$ = uj e - u k - u ab have 
to be large compared to the couplings fl a , £l b , A a , A b , g a 
and g h i.e. |A a |, |A 6 |, |<J*|, |<5 b fe | » |fi a |, |fi 6 |, \g a \, \g b \ and 
|A a |,|Ab| 3> |A a |,|Af,| (for all k), whereas now Raman 
transitions between levels a and b should be suppressed. 
Hence parameters must be such that the dominant 2- 
photon processes are those that involve one laser photon 
and one cavity photon each but where the atom does no 
transition between levels a and b (c.f. figure [5]). To avoid 
excitations of real photons in these processes, we thus 



require |A a - 5%\ , A b - 6%\ > 



2A a 



2A b 



(for all k). 



Whenever two atoms exchange a virtual photon in this 
scheme, none of them does a transition between \a) and 
\b). Moreover both atoms experience a Stark shift that 
depends on the state of the partner atom. This con- 
ditional Stark shifts play the role of an effective a z a z - 
interaction. Dropping irrelevant constants, the resulting 
effective Hamiltonian reads: 



H V r, 



N 

E 



(2) 



where the parameters B and J z are given to second order 
by Here again, the interaction J z and the field B 

can be tuned independently, either by varying Q a and Q b 
for J z or by varying A Q and Af, for B. In particular, |A a | 2 
and |A(,| 2 can for all values of f2 a and Vi b be chosen such 
that cither J z « B or J 2 » B. 

The complete effective model: Making use of the 
Suzuki- Trotter formula, the two Hamiltonians (TTJ) and 
([2]) can now be combined to one effective Hamiltonian. 
To this end, the lasers that generate the Hamiltonian (JTJ) 
are turned on for a short time interval dt ( 1 1 H xy \ \ ■ dt <C 1 ) 
followed by another time interval dt ( 1 1 H zz \\-dt <g; 1 ) with 
the lasers that generate the Hamiltonian © turned on. 
This sequence is repeated until the total time range to 
be simulated is covered. The effective Hamiltonian sim- 
ulated by this procedure is H spin = H xy + H zz or 



N 



H 



spin 



(3) 



a— x,y,z 



where Btot = B + B. The time interval dt should 
thereby be chosen such that 1 , (7 1 <C dt\,dt2 



J- 1 , J- 1 , J~ X ,B- 1 and B' 1 , so that the Trotter se- 
quence concatenates the effective Hamiltonians Hxy and 
Hzz- The procedure can be generalized to higher order 
Trotter formulae or by turning on the sets of lasers for 
time intervals of different length. 

Numerical tests: To confirm the validity of our ap- 
proximations, we numerically simulate the dynamics gen- 
erated by the full Hamiltonian H and compare it to the 
dynamics generated by the effective model ((3]). 

As an example we consider two atoms in two cavities, 
initially in the state -^(l^i) + |&i)) <8 \ a 2), and calculate 
the occupation probability p{af) of the state |oi) which 
corresponds to the probability of spin 1 to point down, 
p{ii)- Figure [3h shows p(ai) and p(li) for an effective 
Hamiltonian © with B tot = 0.135MHz, J x = 0.065MHz, 
J y = 0.007MHz and J z = 0.004MHz and hence \B tot \ > 
\J x \i whereas figure shows p(a\) and p{li) for an ef- 
fective Hamiltonian ([3]) with i3 to t = — 0.025MHz, J x = 
0.065MHz, J y = 0.007MHz and J z = 0.004MHz and 



hence \B t , 
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FIG. 3: The occupation probability p(ai) of state |ai) (solid 
line) and the probability p(|i) of spin 1 to point down (dashed 
line) for the parameters u e — 10 6 GHz, Lo ab = 30GHz, A a = 
30GHz, A,, = 60GHz, tu c = co e - A 6 + 2GHz, A a = 15GHz, 
Q a = Q b = 2GHz, A a = A b = 0.71GHz, g a = g b = 1GHz, 
Jc = 0.2GHz and 5i = -0.0165GHz (plot a) respectively 
Si — —0. 0168GHz (plot b). Both, the occupation of the ex- 
cited atomic states (\ej) (ej \) and the photon number (a" 1 a) 
are always smaller than 0.03. 

Discrepancies between numerical results for the full 
and the effective model are due to higher order terms for 
the parameters 13, 12 1 , which lead to relative corrections 
of up to 10% in the considered cases. Let us stress here 
that despite this lack of accuracy of the approximations 
Tot 12 1 . the effective model is indeed a spin-1/2 Hamilto- 
nian as occupations of excited atomic and photon states 
are negligible. 

Cluster state generation: The Hamiltonian (|2|) can be 
used to generate cluster states [3, [H|. To this end, all 
atoms arc initialized in the states {\aj) + \bj))/y/2, which 
can be done via a STIRAP process [15( , and then evolved 
under the Hamiltonian ([2]) for t = tt/4J z . Figure [4] shows 
the von Neumann entropy of the reduced density matrix 
of one effective spin E v n and the purity of the reduced 
density matrix of the effective spin chain P s for a full 
two cavity model. Since -E v n ~ ^ot^ for t m 19/is while 
the state of the effectve spin model remains highly pure 
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(P s = tr[p 2 ] > 0.95) the degree of entanglement will be 
very close to maximal, see e.g. [lj§. Thus the levels |oj) 
and \bj) have indeed been driven into a state which is, 
up to local unitary rotations, very close to a two-qubit 
cluster state a.k.a. singlet states. 




2:, 



10 15 20 25 5 10 15 

tinl0~ 6 s tinl0~ 6 s 

FIG. 4: a The von Neumann entropy _B v n of the reduced 
density matrix of 1 effective spin in multiples of In 2 and b 
the purity of the reduced state of the effective spin model for 
2 cavities where J z = 0.042MHz. 



evanescent fields or optical fibers that transfer photons 
from one cavity to another. 

Summary: We have shown that single atoms in in- 
teracting cavities that are operated in a strong coupling 
regime can form a Heisenberg spin-1/2 Hamiltonian. All 
parameters of the effective Hamiltonian can be tuned in- 
dividually, making the device a universal simulator for 
this model. When operated in a two dimensional array 
of cavities the device is thus able to simulate spin lattices 
which are not trackable with numerics on classical com- 
puters. Furthermore, this system can be used to generate 
cluster states on such lattices. Together with the possi- 
bility to measure individual lattice sites it thus provides 
the two key requirements for one way quantum computa- 
tion. This demonstrates the versatility of the present set- 
up for the control and manipulation of quantum systems 
in parameter ranges that are experimentally accessible. 



Experimental implementation: For an experimental 
implementation, the parameters of the effective Hamil- 
tonian, J Xl J y , J z , B and B have to bemuch larger than 
rates for decay mechanisms via the photons or the excited 
states \ej}. 

With the definitions = max(f2 a , fij), g = 

max(g a , (?{,), A = min(A a ,A&), the occupation of the 
excited levels \ej) can be estimated to be (\ej)(ej\) sa 
|Q/2A| 2 , whereas the photon number is n p « 
|(£lg/2A)7i| 2 and the couplings J x , J y and J z are ap- 
proximately |(Jlg/2A)| 2 72. Spontaneous emission from 
the levels \ej) at a rate T E and cavity decay of photons 
at a rate Tc thus lead to decay rates Ti = |£7/2A| 2 r£ 
and T 2 = |(%/2A)7i| 2 r c for the effective model. 

Hence, we require Ti <C |(f2g/2A)| 2 7 2 and T 2 <C 
|(rk//2A)| 2 72 which implies r_e < J c g 2 /S 2 and T c < 
Jc (Jc < 5/2), where, S = \{uj a + w&)/2 - uc\ for the 
XX and YY interactions and S = \co — loc\ for the ZZ 
interactions and we have approximated I71I (5 _1 and 
I72I ~ Jc'S~ 2 - Since photons should be more likely to 
tunnel to the next cavity than decay into free space, 
T c <C Jc should hold in most cases. For T E <C J c g 2 /S 2 , 
to hold, cavities with a high ratio g/T E are favorable. 
Since S > 2Jc, the two requirements together imply that 
the cavities should have a high cooperativity factor. 

This regime can be achieved in micro-cavities, which 
have a small volume and thus a high g. Suitable can- 
didates for the present proposal are for example pho- 
tonic band gap cavities [H which can cither couple to 
atoms or quantum dots. Here, cooperativity factors of 
g 2 /2Tc^E ~ 10 and values of g/T E ~ 100 have been re- 
alized and g 2 /2TcT E ~ 10 5 respectively g/T E ~ 10 5 are 
predicted to be achievable [l7j . Further promising de- 
vices are micro-cavities on a gold coated silicon chip that 
couple to single trapped atoms, where g 2 /2TcT E ~ 40 
and g/T E ~ 50 have been achieved [f|. Both are fabri- 
cated in large arrays and couple via the overlap of their 
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